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The CoGeNT collaboration has recently made available new data collected over a 
period of 15 months. In addition to more accurately measuring the spectrum of 
nuclear recoil candidate events they have announced evidence for an annual modu- 
lation signal. We examine the implications of these new results within the context 
of mirror/hidden sector dark matter models. We find that the new CoGeNT data 
can be explained within this framework with parameter space consistent with the 
DAMA annual modulation signal, and the null results of the other experiments. 
We also point out that the CoGeNT spectrum at low energies is observed to obey 
dR/dEji oc which suggests that dark matter interacts via Rutherford scatter- 

ing rather than the more commonly assumed contact (four-fermion) interaction. 
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The CoGeNT experiment operating in the Soudan Underground Laboratory has 
been searching for hght dark matter interactions with a low energy threshold P-type 
Point Contact germanium detector [1, 2]. They have obtained a low energy spectrum 
which is not readily explainable in terms of known background sources, and is con- 
sistent with elastic scattering of light dark matter particles[3, 4, 5]. Recently, 15 
months of data has been analyzed [2] which greatly improves the measurement of 
the low energy spectrum and appears to be annually modulated with a phase con- 
sistent with dark matter expectations [6]. This further strengthens the dark matter 
interpretation of the CoGeNT low energy spectrum. 

The recent CoGeNT results also reinforce the long standing observation of the 
DAMA collaboration of an annual modulation signal in their Nal detector [7, 8]. The 
DAMA signal is extremely impressive with statistical significance of 8.9 sigma with 
phase and period in agreement with the dark matter expectations to high accuracy. 
Attempts to interpret the DAMA signal in terms of a hypothetical background have 
become more and more implausible, and there is ample reason to be confident that 
DAMA and now CoGeNT have observed dark matter. 

The positive results of DAMA and CoGeNT together with the null results of very 
sensitive, but higher threshold experiments such as CDMS[9] and XENON100[10] 
suggest that dark matter is light (~ 30 GeV). Having dark matter light, has long 
been known to alleviate the tension between DAMA and the higher threshold 
experiments[ll] (see also ref.[12]). However, the sensitivity of the higher thresh- 
old experiments has got to the point where there is now some tension between 
DAMA/CoGeNT and e.g. XENONIO. XENONIOO, CDMS/Si, CDMS/Ge when 
interpreted in terms of standard WlMPs even if they are light [13]. Another issue is 
that the allowed parameter regions for the DAMA and CoGeNT signals, although 
close together do not significantly overlap in the standard WIMP framework [5]. It 
turns out that both of these difficulties can be resolved if dark matter is not only 
light, but assumed to be multi-component and self interacting. 

A generic example of this [14] is dark matter from a hidden sector, which contains 
an unbroken U{1)' gauge interaction which is mixed with the standard U{\)y via 
renormalizable kinetic mixing interaction: [15] 



where F^^, is the standard U{1)y gauge boson field strength tensor, and F^^, is the 
field strength tensor for the hidden sector U{1)' . This interaction enables hidden 
sector U{1)' charged particles (of charge Qe) to couple to ordinary photons with 
electric charge Qe'e = ee. We consider the case where the hidden sector contains 
two (or more) stable U{1)' charged dark matter particles, Fi and F2 with masses mi 
and 7712 [Fi and F2 can be fermionic or bosonic]. Under the standard assumptions of 
a dark halo forming an isothermal sphere the condition of hydrostatic equilibrium 
relates the temperature of the particles to the galactic rotational velocity, Vrat' 
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is the mean mass of the particles in the galactic halo. We 
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have assumed that the self interactions mediated by the unbroken 1/(1)' gauge in- 
teractions are sufficiently strong so that they thermalize the hidden sector particles, 
Fi and F2. The interaction length is typically much less than a parsec[16] and 
the dark matter particles form a pressure-supported halo. The dark matter parti- 
cles are then described by a MaxweUian distribution with fi{v) = exp{—E/T) — 
exp{—\miv'^/T) — expl—v^ / VQ{i)\ where 

/ TTX 

voii) = Vrot\ — • (3) 

V m 

With the assumptions that 7712 » rrii and that the abundance of F2 is much 
less than Fi, we have that Vq{F2) <^ v'j^^. The narrow velocity dispersion (re- 
call cr^ = 3^0/2) can greatly reduce the rate of dark matter interactions in higher 
threshold experiments such as XENONIOO whilst still explaining the signals in the 
lower threshold DAMA and CoGeNT experiments. 

While generic hidden sector models are interesting in their own right and have 
been studied in some detail (see e.g. ref.[17]), 1 consider mirror dark matter as 
the best motivated example of such a multi-component self-interacting theory. Re- 
call, mirror dark matter posits that the inferred dark matter in the Universe arises 
from a hidden sector which is an exact copy of the standard model sector[18] (for 
a review and more complete list of references see ref.[19])^. That is, a spectrum of 
dark matter particles of known masses are predicted: e',H',He',0',Fe', ... (with 
rrie' = me,mH' = ttih, etc). The galactic halo is then presumed to be composed 
predominately of a spherically distributed self interacting mirror particle plasma 
comprising these particles [21]. Kinetic mixing of the U{1)y and its mirror counter- 
part allows ordinary and mirror particles to interact with each other and can thereby 
explain the direct detection experiments [3, 4, 11]. The simplest scenario involves 
kinetic mixing induced elastic (Rutherford) scattering of the dominant mirror metal 
component. A', off target nuclei. [The He' and H' components are too light to give 
a signal above the DAMA/CoGeNT energy threshold]. Previous work [3, 4] (see also 
ref.[ll, 14]) has shown that such elastic scattering can explain the normalization 
and energy dependence of the DAMA annual modulation amplitude and also the 
initial (56 day) CoGeNT spectrum consistently with the null results of the other 
experiments, and yields a measurement of e^/^ and rriA'- 



e^U' ~ (7 ± 3) X 10 
rriA' 



-10 



rup 



22 ± 8 (4) 



where ^a' = ua'Ttia' /{0.3 GeV/cm ) is the halo mass fraction of the species A' and 
nip is the proton mass. The measured value of mA'/iTip is consistent with A' ~ O', 
which by analogy with the ordinary matter sector would be the naive expectation. 
Taking a range for ^^/, 1 ~ ^a' ~ 10^^, suggests that e could reahstically range 
from 10"'° to 10"^. Kinetic mixing in this range is consistent with laboratory and 



^Note that successful big bang nucleosynthesis and successful large scale structure requires 
effectively asymmetric initial conditions in the early Universe, T' <^T and rib' /rib ~ 5. See ref.[20] 
for further discussions. 
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astrophysical constraints[22] and has a number of fascinating applications [22, 19, 23]. 
Early Universe cosmology, though, prefers[24] e ~ 10~^. 

The interaction rate in experiments depends on the halo distribution function 
and the interaction cross-section. The former is expected to be a Maxwellian distri- 
bution, fi{v) = exp[—v^/vQ{i)], with vo{i) depending on rfi, as discussed above. In 
the mirror dark matter case, m is expected to be around 1 GeV, but with significant 
uncertainties [3, 4]. Generally it has been found[3] that the dark matter detection 
experiments are relatively insensitive to the precise value of Vq{A') (and hence m) so 
long as Vq{A') <^ v^^^. Kinetic mixing induced elastic Rutherford scattering is par- 
ticularly natural in mirror/hidden sector models as it arises from the renormalizable 
interaction, Eq.(l). In the present study we again assume that it is the dominant 
interaction mechanism coupling ordinary and dark matter. The cross-section for 
a dark matter particle of charge ee to elastically scatter off an ordinary nucleus 
(presumed at rest with mass and atomic numbers A, Z) is given by[ll]:^ 

^'^ ^ (5) 



where 



27re2z2a2 

A = FliqvA) 6 

niA 

and FA{qTA) is the form factor which takes into account the finite size of the nuclei. 
In the case where dark matter particles also have finite size, as in the mirror dark 
matter case, a form factor for those particles also needs to be included. [For elastic 
scattering of mirror nuclei. A' , of atomic number Z' we must replace e — > Z'e in 
the above cross-section formula]. A simple analytic expression for the form factor, 
which we adopt in our numerical work, is the one proposed by Helm[25, 26]. 
The event rate is given by: 



dR 



^ NtUa' -t^ \v\d V (7) 

dEu J\v\>vmi„ dER k 

where Nt is the number of target nuclei per kg of detector and ua' — Pdm^A' /fnA' 
is the number density of halo dark matter particles A' at the Earth's location (we 
take pdm = 0.3 GeV /cm?). Here v is the velocity of the halo particles relative to 
the Earth and v^; is the velocity of the Earth relative to the galactic halo. The 
integration limit, Vmim is given by the kinematic relation: 



{niA + ruA'YER 
'\ 2mAm\, 

The halo distribution function in the reference frame of the Earth is given by, 
fA'{V:-VE)/k = {7rv^[A'])-^/'^exp{-{v+VEy/vl[A']). The integral, Eq.(7), can easily 
be evaluated in terms of error functions[14, 26] and numerically solved. 



^We employ natural units where h = c=l. 
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To compare with the measured event rate, we must include detector resolution 
effects and overall detection efficiency (when the latter is not already included in 
the experimental results): 



where E'^ is the measured energy and ares describes the resolution. The measured 
energy is typically in keVee units (ionization/scintillation energy). For nuclear recoils 
in the absence of any channeling. keVee = keV/q. where g < 1 is the relevant 
quenching factor. Channeled events, where target atoms travel down crystal axis 
and planes, have g ~ 1. In hght of recent theoretical studies [2 7], we assume that 
the channeling fraction is negligible. It is of course still possible that channeling 
could play some role, which could modify the favoured regions of parameter space 
somewhat. 

For this study we consider two of the simplest examples of multi-component 
dark matter models. Following our earlier works[3, 4, 14, 11] we consider mirror 
dark matter with a dominant mirror metal component. A', of atomic number Z'. 
In this case the electric charge of the dark matter particle is eZ'e. ^ The quantity 
Vo{A') is obtained from Eq.(3) with m ~ 1.1 GcV, which corresponds to a He' 
dominated halo, Yhc' — 0.9, expected[28] for e ~ 10~^. We also consider the more 
generic two component Fi, F2 hidden sector dark matter model discussed above, in 
which case vo is less constrained. 

One can define a quantity and compare these theories with experiment. We 
consider the reference point Vrot = 240 km/s which is representative of recent mea- 
surements for the local rotational velocity[29] . The data we consider consists of 
(a) the CoGeNT energy spectrum: 31 bins of width AE ~ 0.05 keVee given in 
the inset of figure 1 of ref.[2]. This spectrum has already been corrected for effi- 
ciency and stripped of background components, (b) The DAMA annual modulation 
energy spectrum in the energy range 2 < E{keVee) < 8. We have taken into 
account systematic uncertainties in energy scale by minimizing over a 20% vari- 
ation in quenching factors, i.e. qnu = 0.30 ± 0.06, qi = 0.09 ± 0.02 for DAMA 
and qce — 0.21 ± 0.04 for CoGeNT. The mirror dark matter candidate provides an 
excellent fit to the data, with X^in/d.o.f values of 23.1/29 for data set (a), and 
8.9/10 for data set (b).^ Favoured regions in the e\/^A')^A' plane can be obtained 
by evaluating contours corresponding to X^{.^\fiA' i^A') = Xmin + 9 (roughly 99% 
C.L. allowed region). In figure 1 we show the parameter regions favoured by the 
data for the Vrot — 240 km/s reference point. The favoured regions for the DAMA 
and CoGeNT signals are in as good an agreement as one might expect given the 
systematic uncertainties \a rneli we have not considered including the fiducial bulk 

"^In our numerical work we allow A' , Z', to have non-integer values, with Z' = A' /2. Since the 
realistic case will involve a spectrum of elements, the effective mass can be non-integer. 

^The low CoGeNT threshold of 0.45 keVee potentially makes the experiment sensitive to the 
e' component via e'-electron scattering, which would be expected to lead to a large rise in event 
rate at low energies [30]. The data is adequately fit by A' — Ge elastic scattering, with no evidence 
for an extra e' — e scattering contribution. This suggests that the e' halo component has a lower 
temperature than the mirror nuclei component. Such a scenario is possible due to the inefficient 
energy transfer between the light e' and much heavier mirror nuclei. 
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volume uncertainty in CoGeNT of ~ 10% and variation of Vrot within its estimated 
~ 10% uncertainty. 




10 20 30 40 50 

Figure 1: CoGcNT and DAMA 99% C.L. favoured parameter (m^/, e\/^A') regions for the 
mirror dark matter model. The reference point Vrot = 240 km/s is assumed. Also shown 
are the 95% exclusion curves evaluated from the null results of XENONIOO, CDMS/Si 
and CDMS/Ge. 

Also displayed in figure 1 is the 95% exclusion limits evaluated for the CDMS/Si [31], 
CDMS/Ge[9] and XENON100[10] experiments^. In computing these limits, we have 
conservatively taken the energy thresholds of these experiments to be 20% higher 
than the advertised values, to allow for systematic uncertainties in energy calibration 
and quenching factor ^. We also show in figure 2a,b, the predicted results for each 
data set for a particular parameter point near the global best fit, as well as a point 
near the best fit for each data set considered separately. 

^There are also lower threshold analysis by XENON10[32] and CDMS[33] collaborations. How- 
ever when systematic uncertainties are properly incorporated, neither analysis is capable of ex- 
cluding light dark matter explanations of the DAMA/CoGeNT signal [34]. 

^Within the mirror dark matter framework the higher threshold experiments such as CDMS/Ge 
and XENONIOO have an important role in probing the heavier ~ Fe' component [35]. 
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Figure 2a: Mirror dark matter versus the CoGeNT spectrum. The soHd hue is for a point 
near the CoGeNT best fit [m^//mp = 26, eV^A' = 5.2 x 10^"*^*^] while the dashed hne is 
for a point near the global best fit [m^'/mp = 24, eV^A' = 5.7 x 10~^^]. 
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Figure 2b: Mirror dark matter versus the DAMA annual modulation spectrum. The 
solid line is for a point near the DAMA best fit [niA' /fnp = 24, &\/^a' = 6.3 x 10""^''] 
while the dashed line is for the point near the global best fit considered in figure 2a 
[niA'/mp = 24, e^/l^ = 5.7 x 10"^°]. 
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It is interesting to compare the 15 month CoGeNT favoured region, as shown 
in figure 1, with results for the same model obtained with the initial 56 days of 
data[3, 4]. The current favoured region is significantly reduced in size. CoGeNT 
data now feature an upper limit on m^' ~ 30 GeV, which is also supported by the 
nuU results of XENONIOO and CDMS/Si. 

The GoGeNT collaboration report [2] evidence for an annual modulation signal 
in their data at about 2.8a G.L. The amplitude of the modulation, averaged over 
0.5 < E{keVee) < 3.0, is measured to be roughly A ^ 0.46 ± 0.17 cpd/kg/keVee. 
This assumes the amplitude and phase are set to theoretical expectations, while 
a larger amplitude is preferred if the phase is left free. For the theories offered 
here, we find that the GoGeNT annual modulation amplitude (averaged over 0.5 < 
E{keVee) < 3.0) is typically around ~ 0.12 cpd/kg/keVee for the parameter region 
near the global best fit, and does not get above 0.20 for any parameter point in 
the global 99% C.L. favoured parameter region (for the reference point Vrot = 240 
km/s). Thus we find an annual modulation somewhat below the CoGeNT central 
value. This difference, though, is not currently, statistically significant, representing 
only a 1.5-2 sigma downward fluctuation from the central value measured in the 
15 X 0.33 month-kg data sample. Obviously future data, especially the measurement 
of the energy spectrum of the annual modulation amplitude, will be important tests 
of the theories considered here. 

Similar results hold for the more generic two component Fi, F2 hidden sector 
dark matter model discussed earlier. For definiteness we have assumed the same vq 
value for F2 (i.e. same rh value) as for A' in mirror dark matter. We have computed 
the before, minimizing over systematic uncertainties in quenching factor. The 

best fit features X^in/d-o.f vahies of 23.0/29 for the CoGeNT data set (a), and 
9.2/10 for the DAMA data set (b). The parameter range favoured by the GoGeNT 
and DAMA data sets (a) and (b) discussed above is given in figure 3 for this case. 
Note that since the electric charge of F2 is ee rather than Z'ee the allowed region is 
shifted c.f. the mirror matter case: e -H- e/Z' . We have also found that the model 
can fit the data for a wide range of VQ{i) values: v^ii) ~ 140 km/s. 

The explanation of the DAMA/LIBRA and GoGeNT experiments considered 
here has a number of interesting features. As noted previously[4] the signals seen in 
these experiments arise predominately from dark matter particle interactions in the 
body of their Maxwellian velocity distribution rather than the tail (as in the model 
of ref.[12, 5]). Because of this, we do not have a great deal of freedom in modifying 
the predicted shape of the spectrum, and thus the agreement of the model with the 
spectrum observed by GoGeNT is a non-trivial test of the theory. In fact in the 
Vq{A') / v'f.^^ — )■ limit, the energy dependence of dR/dEn [Eq.(7)] follows exactly 
that of da/dEpt and is proportional to l/E]^ for Vmin < Vrot and da/dEpt = for 
Vmin > Vrot- [Excepting here the energy dependence of the form factor which is 
relatively minor for Eji ^ 1 ke Vee in germanium] . The 1 / Ej^ dependence of da / dE^ 
follows directly from the masslessness of the exchanged photon in the Feynman 
diagram describing the interaction and is thus a distinctive feature of dark matter 
interacting via Rutherford scattering. For finite VQ{A')/vrot the l/.E'^ behaviour is 



7 



expected provided that Er is sufficiently small that Vmin ~ Vrot, i-e. for 

{mA + mA>y 

For A = Ge, Vrot = 240 km/s and tjia' ~ 18 GeV (the latter suggested by the fit to 
the DAMA annual modulation signal), we have 

dR ^ 1 ^^^^ 
dER E\ 

for Er ~ 1 keVce. This prediction is impressively consistent with the observations 
as indicated in figure 4. CoGeNT's spectrum falls off more rapidly than l/i?^ at 
Er ~ 1 keVee. This suggests the onset of the kinematic threshold Vmin ~ Vrot at 
these energies and is the origin of the m ~ 30 GeV upper limit indicated in figures 
1,3. 
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Figure 3: CoGeNT and DAMA 99% C.L. favoured parameter (m^j) ^a/Cf^) regions for 
the generic hidden sector dark matter model discussed in the text. The reference point 
Vrot = 240 km/s is assumed. Also shown are the 95% exclusion curves evaluated from the 
nuh results of the XENONIOO, CDMS/Si and CDMS/Ge experiments. 



The dark matter explanation of CoGeNT's spectrum offered here can be com- 
pared with the model of ref.[12, 5] which features WIMPS elastically scattering via 
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a contact (four fermion) interaction rather than via Rutherford scattering. The con- 
tact interaction produces a fiat (in Eji) cross-section, excepting the mild (at these 
energies) recoil energy dependence of the form factor. A rapidly falling spectrum 
would then only be expected if dark matter particles are so light that only par- 
ticles in the tail of the halo velocity distribution can lead to recoils with enough 
energy to be observed. In such a scenario the shape of the dR/dEji spectrum nec- 
essarily depends very sensitively on niwimp- Only for nimmp — 7.0 GeV (and with 
standard assumptions) [2] can that model account for the observed spectrum energy 
dependence dR/dEn ~ ^/E]^ at low Er. However the energy dependence is accom- 
modated, rather than explained, which is in contrast to the Rutherford scattering 
scenarios considered here. 




1 I ■ ■ ■ ■ ■ 1 

0.4 ' 
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Figure 4: Low energy CoGeNT spectrum compared with dR/dEn oc where n = 1 

(dashed line), n = 2 (solid line) and n = 3 (dotted line). The data clearly favour the n = 2 
case, which is expected in the mirror dark matter/hidden sector models considered here, 
and is characteristic of dark matter interacting via elastic Rutherford scattering. 



Future experiments should be able to more clearly distinguish mirror/hidden 
sector models from other theoretical explanations - such as the one discussed in 
ref.[12, 5] and many others considered in recent literature - by e.g. precise measure- 
ments of the annual modulation energy spectrum. The mirror/hidden sector models 
predict a characteristic change in sign of the annual modulation amplitude at low 
energies (see figure 6 of rcf. [4]). Distinguishing the mirror dark matter case from 
the more generic hidden sector model might prove more challenging. Whilst the two 
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theories give essentially identical results for the DAMA/CoGeNT experiments mir- 
ror dark matter predicts a spectrum of particles of known masses. In particular the 
scattering of low mass components, e' on electrons and He' / H' on target nuclei can 
ultimately be seen in very low threshold experiments. Higher mass sub-components, 
such as a ~ Fe' or ~ Ca' component would also be expected and should ultimately 
be observed if dark matter is of the mirror type. 

In conclusion, we have examined mirror/hidden sector dark matter in the light 
of CoGeNT's more precisely measured spectrum and annual modulation signal [2]. 
The CoGeNT spectrum is observed to obey dR/dEn oc ^/E\ at low energies which 
suggests that dark matter interacts via a massless or light mediator (Rutherford 
scattering) rather than the more commonly assumed contact (four-fermion) inter- 
action. Such Rutherford scattering is a feature of mirror and more generic hidden 
sector dark matter models considered here and in previous works[3, 4, 14, 11]. We 
have found that such models provide an excellent fit to the data which is easily 
consistent with the null results of the sensitive but higher threshold experiments, 
such as CDMS and XENONIOO. The next generation P-type Point Contact detec- 
tors, including CoGeNT(C-4), MAJORANA, GERDA and CDEX should be able 
to provide a decisive test of these models by e.g. a precise measurement of the an- 
nual modulation energy spectrum. These and other experiments are awaited with 
interest. 
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